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ON A MODEL BY MUMFORD AND NITZBERG 
FOR IMAGE SEGMENTATION 


G. BELLETTINI(') AnD M. PAOLINI(*) 


Abstract. We study a model of image segmentation recently suggested by Mumford and Nitzberg 
which involves the length and the curvature of the contours . The model is described by an energy 
functional defined on ordered families of subsets of R?. More precisely, let g E L°(R?) be a given 
function with compact support, representing the grey-level of the image. The aim is to reconstruct g by a 
minimum point of the functional P defined by 

n n n+1 
P(E1,...,En)= 134) [1+ V(0;,1i)] dH +Dia+d f lgi — gl? dz. 
i=1 S2Ei ii VPI 


s=1 i=l 


Here n € N,n> 1, E1,...,En € R? are bounded open sets of class C?, 9;(z) is the angle that the tangent 
unit vector of dE; at 2 € 9E; forms with the positive direction of the r-axis, x;(z) is the curvature of 2E; 
at 2, |E;| is the Lebesgue measure of E;, H! denotes the one-dimensional Hausdorff measure in R?, and 


i-1 n 
Ej=E, Ei=E\|]E vi=2...,n, E = ]zn 
j=1 g=1 
1 ; 
Mise gdz Va=1;- 9% In+1= 0. 
E Ja 
Finally w: Rx R + [0, +00] is a Borel function satisfying suitable convexity and growth conditions. The 
fact that the functional is defined on partitions of sets and the presence of a (possibly linear) curvature 
term in the energy functional permit to reconstruct the hidden parts of the image occluded by regions 
which are closer to the observer. To deal with corners and cusps, which are viewed as points where the 
curvature is concentrated, we introduce the class of the continuous curves having derivative of class BV, 
and we study their arguments. The proof is based on the techniques of relaxation theory in the Calculus 
of Variations. We estimate the lower semicontinuous envelope of the energy functional, and we compute 
it in some special cases. 
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SU UN MODELLO DI MUMFORD E NITZBERG 
PER LA SEGMENTAZIONE DI IMMAGINI 


G. BELLETTINI & M. PAOLINI 


1. Introduzione. In questo lavoro illustriamo alcuni risultati ottenuti in [4] relati 
vamente a un modello proposto da Mumford e Nitzberg sul problema della segmentazione 
di immagini [26,29]. L’approccio è variazionale, e consiste essenzialmente nello studio 
della semicontinuità di un funzionale dipendente dalla lunghezza e dalla curvatura del 
bordo di un sottoinsieme regolare limitato di R? (si veda la sezione 3). Per tenere conto 
della profondità relativa degli oggetti, il funzionale che descrive il modello viene definito 
su famiglie ordinate finite di sottoinsiemi di R? di classe C2 (si veda (3.1)), e poi esteso 
per inferiore semicontinuità su famiglie formate da insiemi anche non regolari. Il fatto di 
definire un funzionale su famiglie ordinate di insiemi è una delle due novità dell’approccio 
di Mumford e Nitzberg: invece di basare il modello su un insieme di curve che dividono 
un certo dominio in regioni disgiunte [27], si cerca un insieme di porzioni del piano la cui 
unione sia il dominio, e che possano avere zone di sovrapposizione, porzioni legate tra loro 
dall’ordinamento dato dalla profondità. Le seconda novità consiste nel fatto che il termine 
di curvatura ha crescita lineare all’infinito. Ciò distingue in modo essenziale questo mod- 
ello da quello studiato in [3], in cui viene considerato un funzionale (sempre in dimensione 
due) con un termine di curvatura a crescita p > 1 (si vedano anche i lavori [2,20,21,25]). Il 
fatto che il termine di curvatura possa essere a crescita lineare porta a considerare famiglie 
di curve chiuse di classe W1:1 (parametrizzate con velocità costante) aventi come derivata 
una funzione di classe BV (Definizione 2.4). Su tali curve, o più precisamente, sulla classe 
delle funzioni angolo che il vettore tangente a tali curve forma con la direzione positiva 
dell’asse delle x (si veda Lemma 2.1), si può definire un funzionale (si veda (3.7)) che sem- 
bra essere un buon candidato per lo studio del rilassamento della parte di curvatura del 
funzionale originale, che invece era definito su insiemi (si veda (3.5)). Definire il funzionale 
su angoli comporta alcune difficoltà, viste le varie identificazioni modulo 2r di cui bisogna 
tenere conto (si veda la definizione dello spazio BV(]a,b(;R/27)). D'altra parte, come si 
vedrà nel seguito, pensare i funzionali dell’energia definiti su angoli sembra essere la cosa 
più naturale per valutare energie concentrate su punti singolari, come angoli o cuspidi (si 
veda anche la diseguaglianza 2.2). Questo corrisponde a definire il funzionale su funzioni 
a valori nella sfera unitaria S! [18]. 
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Quello che interessa qui è lo studio delle proprietà di semicontinuità e rilassamento 
delle energie (si veda (4.1) e sezione 5). In particolare, oltre alla semicontinuità dei fun- 
zionali sugli insiemi regolari di classe C?, ci interessa valutare l’inviluppo semicontinuo (o 
funzionale rilassato) rispetto alla convergenza in misura di insiemi, in modo da potere dare 
un significato alla curvatura nei punti singolari. Come vedremo in seguito, tale funzionale 
tiene conto delle energie concentrate in tali punti, ma anche delle mutue distanze tra essi. 

Per quanto riguarda la bibliografia su problemi della segmentazione affrontati con 
tecniche variazionali, si vedano tra gli altri [1,5,7,11,13,14,17,22,24,27,28,30,31,32]. 


Infine, per ulteriori dettagli e per le dimostrazioni si rimanda a [4]. 


2. Alcune notazioni e risultati preliminari. Denotiamo con H” la misura di 
Hausdorff h-dimensionale in R?, per A = 0, 1 [16]; indichiamo con |B| la misura di Lebesgue 
di un insieme di Borel B C R?. 

Denotiamo con C?(R?) la famiglia di tutti i sottoinsiemi di R? aperti limitati di classe G°. 
Nel seguito indichiamo con / un intervallo aperto e limitato di R. 

Con (7) = r([0,1])) = {(t) : + € [0,1]} denotiamo la traccia di una curva y e con (7) la 
sua lunghezza. Inoltre, se 2 € R? \ (7), I(Y, 2) è l’indice di 2 rispetto a 7 [6], [9, III.8]. 


DEFINIZIONE 2.1. Un sistema di curve è una famiglia finitaT = {y},...,} di curve 
chiuse di classe W11 tale che per ogni î = 1,...,m la quantità |4Î| sia costante quasi 
ovunque su [0, 1). La traccia (I°) del sistema I si definisce come Ut, (YÎ), e la lunghezza 
IT) come VT, IYÎ). Se 2 € R? \ (T), l’indice I(T, 2) di 2 rispetto a I° è definito come 
3: I(9% z). 


Il sistema di curve D = {y},...,7"} si dice disgiunto se (7°) N (77) = 0 per ogni 
ideati 

Sia E € C?(R?); diciamo che il sistema disgiunto di curve I è una parametrizzazione 
orientata di E se ogni curva del sistema è semplice, (I) = E, e 


E={zeR?:I(F,z)=1}, R?\E={zeR?:I(1,z)= 0}. 


2.1. Gli spazi BV(I) e BV(]a,b[;R/27). Indichiamo con BV(I) lo spazio delle funzioni 
a variazione limitata su I. Per le principali definizioni e proprietà delle funzioni BV si 
rimanda a [10,16,19,23,33]. 

Data f € BV(I), scriveremo 

f=efatf, 

dove f° € L!(I) è la densità della parte assolutamente continua di f rispetto alla misura 
di Lebesgue, e f? è la parte singolare. 

Se f € BV(I) indichiamo con Sy l’insieme di salto di f, e poniamo f(t-) = ap— 
liminf f(7) S f(t+) = ap — limsup f(7). 
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Identificheremo talvolta f € BV(I) con il suo rappresentante ft) = (f(t-)+f(t+))/2 
definito ovunque (si veda [16, Th. 4.5.9]). 


Dato che intendiamo lavorare su funzioni angolo, introduciamo la seguente definizione. 


DEFINIZIONE 2.2. Siano 91,0, € BV (]a,b[). Diciamo che 0, e 0, sono equivalenti se 
per quasi ogni t € ]a, b[ esiste k(t) € Z tale che 01(t) = 02(t) + 2k(t)m. Indichiamo con [0] 
la classe di equiavalenza di 0, e poniamo 


BV(Ja,b;R/27) = {[9]:0 € BV(]a,5)}. 


È facile vedere che dato [6] € BV(]a,b[;R/27) esiste una funzione © € [0], che chi- 
ameremo un rappresentante minimale di [0], tale che 


(2.1) 1O/(1a, 60) = min{[6l(a,b]) : 8 € [}. 


Denotiamo con M[0] l’insieme dei rappresentanti minimali © di [9] € BV(]a,b[ ;R/27) 
tali che O(a+) € [0, 27]. 


DEFINIZIONE 2.3. Sia [9] € BV(]a,b[ ;R/27), e sia {[0x]}x € BV(]a,b[;R/27). Di- 
‘ciamo che [9n] + [9] per h + +00 se, per ogni @ € [0] e ogni h € N, esiste 0, € [0] tale 
che 9, + @ in L'(]a,b[) perh + +00. 


2.2. Sistemi di curve con derivata discontinua. 


DEFINIZIONE 2.4. Diciamo che un sistema di curve D = {a1,...,7} è di classe B se 
gi e W14([0, 1]; R?) e ji € BV([0,1]; R?) per ogniî= 1,...,m. 


Sia y una curva di classe 8; se B C [0, 1] è un sottoinsieme di Borel che contiene 
uno degli estremi dell’intervallo [0, 1], con |5|(B) intendiamo I#I(Bn ]0, 1) +|XI({0}), dove 
14I({0}) si considera definito per l’estensione periodica di 7. 


2.3. Limiti di sistemi di curve regolari. 


DEFINIZIONE 2.5. Diciamo che una successione {Ta} di sistemi di curve di classe C2 
è debolmente convergente a una famiglia di curve D = {y},...,y°} se il numero dell curve 
di ogni sistema Tx è lo stesso del numero di curve di T per h sufficientemente grande, i.e., 
Ta = {M--- 97°} per ogni h, e, inoltre, Xi + y' uniformemente, fi — 4* debolmente in 
BV per h + +00, per ognii= 1,...,m. 


x 


Osserviamo che se una successione di sistemi di curve {T1}x di classe C? è debolmente 
convergente alla famiglia D = {y},...,"}, allora D è un sistema di curve di classe B e 
I(1, 2) € {0,1} per ogni 2 € R? \ (1). 
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DEFINIZIONE 2.6. Diciamo che I è un sistema limite di curve se I° è il limite debole 
di una successione di parametrizzazioni orientate di aperti limitati di classe C, 


2.4. Esistenza della funzione argomento per curve con derivata discontinua. Il seguente 
lemma permette di parlare di argomento di 4, dove y è una curva di classe B parametrizzata 
con velocità costante. 


LEMMA 2.1. Sia f € BV(]a,b[ ;R?) tale che |f(t)| = R per quasi ogni t € ]a,b[, per 
un certo R > 0. Allora esiste una funzione © : ]a,b[ + R che ha le seguenti proprietà: 
(a) © € BV(]a,b[) e f(t) = (Rcos0(t), Rsin O(t)) per quasi ogni t € ]a, bl; 

(b) Se= Sf; 
(c) —n < O(t+) — O(t—) < rr per ogni t € Ja, bl. 


Quando f = #4 e 7 € B, ci vuole una attenzione ulteriore dovuta al fatto che la curva è 
chiusa, e dunque parametrizzata su un intervallo con estremi identificati; pertanto bisogna 
specificare cosa si intende per variazione del relativo argomento su tale intervallo. Se 
B C [a,b] è un insieme di Borel che contiene almeno uno degli estremi dell’intervallo di 
parametrizzazione, con |Ò|(B) intendiamo |O|(BN ]a,b[) + |O(a+) — O(6—) — 2k7], dove 
©(a+) = ap — lim ©(#) = lim O(t), O(6-) = ap— lim ©(t) = lim O(t), ek € Zè tale 

t+at t+at t4+b7 tab 
che —r < O(a+) — O(b—) — 2kr < 7. 


Una osservazione utile è la seguente: se © è un argomento di 4, e y € B, si ha 
(2.2) 1521({s}) <IO°I({s}) Vs € Se. 


Questa diseguaglianza mostra che lavorare su curve non equivale a lavorare su funzioni 
angolo, dato che le energie concentrate su punti di discontinuità sono differenti. Questo 
discorso si collega a quanto detto nell’introduzione. Abbiamo scelto le funzioni angolo 
perchè ci sono sembrate più naturali in vista dei risultati che uno vuole ottenere (come la 
penalizzazione dei punti singolari del bordo di un insieme in termini del salto dell’angolo). 

La seguente proposizione mostra il legame tra le convergenze delle curve e quelle dei 
relativi angoli formati dai vettori tangenti. Infatti, supponiamo di avere un insieme E 
regolare tranne che in un punto di cuspide, e consideriamo la successione approssimante di 
insiemi che coincide con E tranne che in un intorno del punto di cuspide, dove è costruita 
con archetti di cerchio sempre più piccoli che si stringono intorno alla cuspide (si confronti 
con (3.6)). Se si considerano i relativi angoli formati dai vettori tangenti alle curve, ci 
possono essere scelte dei rappresentanti nelle classi di equivalenza per cui gli angoli delle 
curve approssimanti non convergono all’angolo originale. Tuttavia si può dimostrare il 


fatto seguente. 


PROPOSIZIONE 2.1. Sia y € B, sia {Yn}® C B una successione di curve tale che n + % 
in L([0, 1]; R?) per a + +00. Supponiamo che Y e ciascuna curva Yh siano parametrizzate 
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con velocità costante su [0, 1]. Siano 0,0% € BV ([0, 1]) argomenti di 4,4 rispettivamente, 
per ogni h € N. Allora [9] + [0] per h + +00. 


2.5. Il funzionale rilassato. Dato un funzionale £ a valori reali estesi e definito sulla 
classe degli insiemi misurabili di R?, indichiamo con £ l’inviluppo semicontinuo (o fun- 
zionale rilassato) di £ rispetto alla topologia di L!(R?) (si veda [8]). 


3. I funzionali che descrivono il modello. Sia 9g € L°(R?) una funzione asseg- 
nata a supporto compatto. Tale funzione rappresenta l’intensità di luce prodotta da una 
fotografia. Lo scopo è quello di riscostruire 9, 0, come si usa dire, segmentare l’immagine g. 
Questo significa trovare regioni di particolare interesse dell’immagine stessa, che possano 
servire per darne una descrizione approssimata. In questo modello tali regioni interessanti 
possono sovrapporsi e avere zone di occlusione e, teoricamente, le parti nascoste possono 
essere in qualche modo riscostruite. 

Vogliamo dunque riprodurre g considerando un punto di minimo del funzionale P 
definito come 


n+1 


(3.1) PEis....B)=Y f_ [1+4(09;,x:)] dA! +D12+Y f, Igi — gl? de. 


i=1 s=1 


In questa formula n € N, n > 1, Ei,...,En € C}(R?), 0i(2) è l'angolo che il versore 
tangente a VE; nel punto 2 € 0E; forma con la direzione positiva dell’asse delle x, Ki(z) è 
la curvatura di 0E; in z, e 


i-1 n 
Ei=E,  Ei=E\\JE; Vi=2,...,n, E+41=R?\{]E;, 
(3.2) j=1 j=1 
1 
9= TE] dia Vie Last In+1 = 0. 


Infine w : R x R + [0, +00] è una funzione di Borel che soddisfa le seguenti ipotesi: 
(i) (7, -) è convessa su R per ogni 7 € R; 

(ii) + è inferiormente semicontinua su R x R; 

(iii) %(-,0) è limitata; 

(iv) esistono due costanti positive A > 0 e B > 0 tali che 


(3.3) v(n,é) > AléÉ][-B V(né)ERxR; 


(v) (n, €) = %(2+ 7, é) per ogni (7,6) €R x R; 
(vi) v(n,é) = W(7,—-é) per ogni (n,é6) ERxR. 
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Le ultime due assunzioni su # servono per semplificare la trattazione, ma alcuni risultati 
resterebbero validi anche senza di esse. Come caso modello per % si può considerare la 
funzione w definita come 


cié? se |E| < T, 
34 sz di 
(3.4) W(n,é) = %(I€|) { cal] — c3 se |É|>T, 


con T > 0, ca > 0, cir? = c2T—c3, c2 = 2c17; questa è la funzione adoperata da Mumford 
e Nitzberg per la descrizione del modello. 

Quello che si deve pensare è che ogni partizione {E,..., En} è munita dell’ordinamento 
seguente: E; > E; (i.e., E; è più vicino di E; all’osservatore) se è < J. Questo ordinamento 
rappresenta la profondità relativa (o l’occlusione); dunque, se è = 1,...,n, l'insieme Ei 
definito in (3.2) è la parte visible di E;, e En+1 è lo sfondo. 

Il termine in ? contenente la funzione g serve per mantenere g pressochè costante su 
ogni regione. Il termine di area penalizza l’area totale delle regioni, e incoraggia segmen- 
tazioni in cui piccole regioni occludono regioni più grandi. Infine, il termine di lunghezza 
e curvatura richiede che le frontiere delle regioni abbiano lunghezza non troppo grande e 
non siano troppo curvate. 

Poichè siamo interessati alla semicontinuità di P rispetto alla convergenza I} di insiemi, 
possiamo semplificare il funzionale trascurando subito due termini (continui rispetto alla 
topologia L'), quello dipendente dalla g e quello dipendente dall’area. Inoltre non è difficile 
vedere che tutto si riduce a considerare il funzionale semplificato 


F(E) = 1) It] e. 


Nel seguito dunque ci concentreremo sul funzionale F. Pertanto precisiamolo meglio. Sia 
E € C2(R?), e sia T = {y},...,y7} una parametrizzazione orientata di GE. Sia 0(2) 
l'angolo che il versore tangente a VE in 2 € 8E forma con la direzione positiva dell’asse 
delle x. e sia #(2) la curvatura di VE in 2. Definiamo 


H!(0E) + i (0,8) dH! se E € C£(R?), 
(3.5) F(E)= ai 


+00 altrove sugli insiemi misurabili. 


Poniamo F(09) = 0. 

Semplici esempi mostrano che esistono insiemi non regolari E tali che F(E) < +00, 
come i poligoni e gli insiemi con cuspidi. È facile rendersi conto che, se uno ha un insieme 
regolare tranne che in un punto p di angolo, si può trovare una successione di insiemi 
regolari che approssima E in area e al limite fornisce una energia finita. Infatti, prendiamo 
come in (3.4), e indîthiamo con a l’angolo (interno) formato nel punto p dalla frontiera 
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di E. Basta allora considerare la successione formata da insiemi che sono uguali ad E 
dappertutto tranne che in un intorno di p, e ivi hanno un archetto di circonferenza che 
smussa il punto in questione, e ivi si stringe. Questo fornisce immediatamente la stima su 
F(E) seguente: 


(3.6) F(E) < H!(0E) + dai i W(|k|) dH' + c2(7 — a). 


Osserviamo che il punto singolare p è dunque penalizzato, in termini dell’energia, di un 
fattore proporzionale a 7 — a. Tale penalizzazione corrisponde alla variazione dell’angolo 
del versore tangente alle curve approssimanti, in un intorno di p. 


Le domanda naturali sono allora le seguenti: come definire un funzionale su insiemi non 
regolari, e come rappresentarlo usando le parametrizzazioni dei bordi? L’estensione, come 
comunemente si fa in Calcolo delle Variazioni, si fa usando F. La scrittura di tale funzionale 
come funzione di insieme risulta invece un problema, a nostra conoscenza, ancora aperto. 


Introduciamo per cominciare i funzionali su funzioni angolo e su curve aventi discon- 
tinuità. Indichiamo con +oo la funzione di recessione di + rispetto alla variabile E, i.e., 


Yoo(n, 6) = lim t4(n,6/t) V(né eRxR. 
t+0+ 
Sia 9 € BV(]a,b[;R/2r), 0 € [9] e © € M[6]; per ogni intervallo aperto / C Ja, b[ definiamo 
(3.7) K(9,I)= K°(9,I1)+K°(0,1), 


dove 


K°(0,I)= J v(4,8°) dt, 
I 


O(s+) 
/ Woo(7, 1) dr. 


G9 EGD=/  Be(01) A+ Li 


sEINSe 


Nella espressione di K* abbiamo O(s-) < 0(s+), dato che O(st) sono i limiti inferiore e 
superiore approssimati di © nel punto s. Si può mostrare che K°(9,I ) non dipende dalla 
scelta di @ € [9], e che K*(9,1) non dipende dalla scelta di © € Md]. 

Se © è un argomento di 4 per una opportuna curva y : [a,b] + R? di classe B 
parametrizzata con velocità costante, poniamo 


K°(0, [a,b]) = K*(0,]a, b[) + 


O(a+)-2kr 
Î Woo(7,1) drl, 
CICHS) 


dove k € Z è tale che —r < O(a+) — O(b-) -— 2kr < rx. 
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Si può dimostrare che il funzionale appena definito risulta semicontinuo rispetto alla 
convergenza L? (delle classi di equivalenza). 


Definiamo adesso il funzionale sulle curve. Sia I = {},...,} un sistema di curve 


di classe B parametrizzato per lunghezza d’arco, sia ©; un argomento di ' su [0, (7°)] per 
ogni i= 1,...,m. Definiamo 


K(F)=K°(1)+K°(T),  F(T)=M)+K(1), 


dove 


K°(1) = YO K"(0;,0,10")]),  K*M=) "(00 {yÎ))). 
i=l 


i=l1 


4. Semicontinuità di F. Vale il seguente teorema, di cui vogliamo commentare la 
dimostrazione. 


TEOREMA 4.1. Sia E € C?(R?), e sia {En}n C?(R?) una successione convergente 
ad E in L'(R?). Allora 


(4.1) F(E) S liminf F(Fx). 


‘TRACCIA DELLA DIMOSTRAZIONE. Possiamo supporre che il membro destro di (4.1) 
sia finito, altrimenti non c’è niente da dimostrare. Sia {Ex}: una sottosuccessione {En} 
tale che 


pi PB) = mint (E) < 400 


Per semplicità, questa successione (e ogni ulteriore sua sottosuccessione) verrà denotata 
con {Ex}:. Grazie al fatto che le energie degli insiemi approssimanti sono equilimitate 
rispetto a h, si può mostrare che, per ogni k, 0Ex ha un numero finito di componenti 
connesse. Sia Ax una parametrizzazione orientata di 0Ex. Sempre per motivi di equi- 
limitatezza delle energie, possiamo supporre, passando a eventuali sottosuccessioni, che 
Ax = fol, .-.-,90E} conm indipendente da £. A questo punto si riesce a dimostrare, con 
argomenti di compattezza, che esiste una successione {Tx}: = Di, sa sr} (il cui ele- 
mento k-esimo è un sistema di curve ottenuto da foi, ...:7°} eventualmente trascurando 
quelle curve con lunghezza infinitesima, o che hanno traccia che si allontana all’infinito 
rispetto all’origine di R?), ed esiste un sistema limite di curve I € B, che verificano le 
seguenti proprietà: n £ m, PS Ty, af e, per ogni j = l,...,n si ha pi + gli 
uniformemente e 4} — 4ii debolmente in BV per k + +00. Denotiamo con 0} (rispet- 
tivamente 0;,) un argomento di sj (rispettivamente di 45). Dalla Proposizione 2.1 segue 
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che oi + @;; in L! per k + +00. Abbiamo dunque 


pm) = dm FC) > de D0Ò) 


n pl) n'e LI 
+ liminf 3 J v(OÙ, 6h) ds = x (95) 
e " 1 $ s DICA 
+ mia Yo) IO RTORCE” 
n 1 - i . 
> i#\limi i; (y5\-16% 
210) + DO 9*) ii frei) a. 
Ma i funzionali sugli angoli si dimostrano essere semicontinui. Ne deduciamo che 
liminf F(E4) > IT) +) K(©:,,[0,(95)]) = {M)+K(0) = F(P). 


a=1 


Adesso si riesce a vedere che vale sempre l’inclusione (T) > dE. Inoltre, vale anche un 
lemma che permette di paragonare le energie di I e E quando (T) D dE. Precisamente, 
si dimostra la cosa che segue: sia E € C$(R?), sia T € B un sistema di curve tale che 
(1) > dE. Allora 


J v(0,) dH! < K°(T). 
dE 
Ne deduciamo F(T) > F(E), e la dimostrazione è conclusa. O 


5. Qualche esempio di calcolo di F(E). In questa sezione la funzione w= %(|€]) 
è quella in (3.4). Si può dimostrare il seguente teorema, che chiarisce la diseguaglianza 
(3.6). 


TEOREMA 5.1. Sia E un insieme regolare tranne che in un punto p, dove dE ha un 
angolo (interno) a, con 0 < a < 2r. Allora 


H'(0E) + J %(|k|) dH! + ca min(a, 27 — a, |r — |) 
dE\{p) 
<F(E)<H'(0E) + fi W(|e]) dl + calm — al. 
9E\{p} 


TRACCIA DELLA DIMOSTRAZIONE. Con un ragionamento come quello fatto in (3.6), ci 
possiamo ridurre a dimostrare che 


F(E)> H! (0E) + J %(lk|) dH! + c2 min(a,2r — a, |T— al). 
9E\{p} 
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Grazie ad opportuni risultati di compattezza, si dimostra che 
FE) > W(0E)+ f__W(lel) dl + FT, {D}), 
dE\{p} 


dove IL = {},...,"} è un sistema limite di curve di classe B la cui traccia contiene dE, 
e tale che E coincide con l’insieme dei punti di indice 1 rispetto a I; inoltre, se ©; è un 
argomento di 4' su [0, I(7°)], abbiamo definito 


re.=L XL f 


i=1 seyi tp)” 94) 


©;(s+) m 
Wo(7,1) dr=c2)} )} |0i(s+)-9i(s-)| 


i=1 ser! (p) 


Questa grandezza descrive l’energia concentrata del sistema I nel punto p. Si tratta dunque 
di dimostrare che 
F(T;{p}) > ca min(a, 27 — a, |r — al). 


Questa diseguaglianza si dimostra in vari passi, usando le proprietà di parità dei vettori 
tangenti a I nel punto p [4]. Tutto è basato sul fatto che E = {: € R? :I(1,2)= 1}, che 
è a sua volta conseguenza del fatto che I è limite di bordi di insiemi. D 


COROLLARIO 5.1. Se alle ipotesi del Teorema 5.1 si aggiunge 7/2 < a £ 37/2, allora 
(5.1) F(E)=H'(0E)+ J W(|x]) dH! + co|r — al. 
9E\{p} 


Noi pensiamo che quest’ultimo risultato valga senza restrizioni su a. Il motivo è che 
non ci sono altri punti di singolarità diversi da p in cui una linea nascosta uscente da p 
possa terminare. 

Sottolineiamo che una eguaglianza tipo la (5.1) dove però siano presenti più punti di 
angolo a; (ognuno che dia il suo contributo di energia concentrata) non vale in generale, 
se a; < T/2 0 a; > 37/2 per qualche i. Vale infatti il seguente teorema, che mostra che 
in certe condizioni l’energia minima possibile (i.e., il valore del funzionale rilassato) non è 
raggiunta smussando gli angoli con archi di cerchio, ma collegando due angoli tra loro con 
dei tubicini sempre più stretti (si confronti anche con [3]). 


TEOREMA 5.2. Supponiamo che E sia l’unione di due componenti connesse E, ed Ep, 
a forma di goccia, affacciate l’una all’altra, aventi ciascuna un solo punto singolare dato 
da uno stesso angolo interno 0 < a < r/2. Siano pi e pz i punti singolari di E e Ea, 
rispettivamente, e supponiamo che tali punti siano allineati, e siano proprio i punti da cui 
le due gocce “si affacciano”. Se la distanza D tra pi € p2 è sufficientemente piccola, allora 


F(E)=#!(0E)+ J, v(|x]) dH! +2D + 2c20 
dE\{p1,p2} 
<H(0E) + J, 4(|x|) dH! +2c2(r — a). 
0E\{p1,p2} 


83 


TRACCIA DELLA DIMOSTRAZIONE. Anche in questo caso, si riesce a dimostrare che 


FE) = #(0E)+ f eos DI 
P1,P2 


#0} [1+6068)] +0), 
i=1 ‘{s€0,Uy)]:r' (3) 60E} 


dove I è un sistema di curve con proprietà simili a quelle elencate nella dimostrazione del 
Teorema 5.1. Resta da dimostrare che, se D è sufficientemente piccolo, allora 


m 
J [1 + v(1621)] ds+K*(T)>2D+2cz0. 
i=1 V{s€[0,Uvi)]:vi (9) €0E} 


Anche in questo caso bisogna usare argomenti di parità sui vettori tangenti a I nei punti 
Pi € pz, un poco più complicati di quelli adoperati per dimostrare il Teorema 5.1. O 

Concludiamo con un breve commento al modello esposto. Tale modello ha alcune 
limitazioni (si veda la discussione in [29, Ch. 6]): in particolare, notiamo che viene imposto 
sostanzialmente un solo valore di intensità di g in ogni regione che forma la segmentazione. 
Tale limitazione può essere superata considerando, al posto del termine i-esimo in cui 
compare g nella scrittura di P, un termine tipo 


J (ou? dr + J Vul? de, 
E: Ei 


dove la funzione u; approssima l’intensità di luce osservata nella regione E; ed è ivi regolare. 
Naturalmente con questa differenza lo studio del modello diventa più complicato. 
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